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ABSTRACT 

The analysis of the rotation curves (RCs) of spiral galaxies provides an efficient diagnostic for 
studying the properties of dark matter halos and their relations with the baryonic material. Since 
the cored pseudo-isothermal (Iso) model usually provides a better description of observed RCs than 
does the cuspy NEW model, there have been concerns that the ACDM primordial density fluctuation 
spectrum may not be the correct one. We have modeled the RCs of galaxies from The HI Nearby 
Galaxy Survey (THINGS) with the Einasto halo model, which has emerged as the best-fitting model 
of the halos arising in dissipationlcss cosmological A^-body simulations. We find that the RCs are 
significantly better fit with the Einasto halo than with either Iso or NEW halo models. In our best-fit 
Einasto models, the radius of density slope —2 and the density at this radius are highly correlated. 
The Einasto index, which controls the overall shape of the density profile, is near unity on average 
for intermediate and low mass halos. This is not in agreement with the predictions from ACDM 
simulations. The indices of the most massive halos are in rough agreement with those cosmological 
simulations and appear correlated with the halo virial mass. We find that a typical Einasto density 
profile declines more strongly in its outermost parts than any of the Iso or NEW models whereas it 
is relatively shallow in its innermost regions. The core nature of those regions of halos thus extends 
the cusp-core controversy found for the NEW model with low surface density galaxies to the Einasto 
halo with more massive galaxies like those of THINGS. The Einasto concentrations decrease as a 
function of halo mass, in agreement with trends seen in numerical simulations. However they are 
generally smaller than values expected for simulated Einasto halos. We thus find that the Einasto 
halo model provides, so far, the best match to the observed RCs, and can therefore be considered as 
a new standard model for dark matter halos. 

Subject headings: cosmology: dark matter - galaxies: halos - galaxies: structure - galaxies: spiral - 
galaxies: kinematics and dynamics - galaxies: fundamental parameter 



I. CONTEXT 

One of the major tools in galactic dynamics is the de- 
composition of the rotation curves (hereafter, RCs) of 
spiral and lenticular galaxies into baryonic and dark com- 
ponents. This mass decomposition produces constraints 
on the distribution of dark matter in these galaxies, 
which permits the investigation of the possible relations 
between baryons and dark matter, and most importantly 
to study the properties of galactic dark matter halos. 

Two spherical halo models are usually used for the 
mass decomposition of RCs. On one hand, cosmologists 
tend to favor the Navarro et al. (1996, hereafter, NEW) 
model 

PNFw(r) = 4p_2 — ( '^'^ ) , (1) 

where r_2 is the radius where the density profile has a 
(logarithmic) slope of —2 (the "isothermal" value) and 
/5_2 is the local density at that radius. The NEW model 
was originally thought ()Navarro. Frenk fc Wliit^ [1991 
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to provide a universal description of halos of different 
mass produced in dissipationlcss (i.e. dark matter only) 
cosmological simulations run with different cosmologies, 
in particular in the Lambda Cold Dark Matter (here- 
after, ACDM) paradigm. It has a cuspy inner structure 
with an inner slope of — 1 and logarithmically diverges in 
mass at large radii. Note that steeper inner cusps (inner 
slopes as steep as —3/2) have been pro duced in other dis- 
sipationlcss cosm(3logical simulations (iMoore et al.|[T999l : 
iJing fc Sutol[2000l:IDiemand et al.ll200l ! 

On the other hand, the analysis of RCs has tradition- 
ally used the so-called pseudo-isothermal (hereafter, Iso) 
model 

a? 

Piso(0 ==4p(a) „ ^ , (2) 

which is often called the core halo because of its finite 
central density (but linearly diverging mass at large radii, 
where the asymptotic density slope is —2). This halo has 
no cosmological background but is often seen to better fit 
galactic RCs than the NEW model, particularly for dark 
matter d ominated objects l i ke low surface density galax- 
ies (e g. Ide Blok fc Bosmal 120021 : iKuzio de Narav et al.l 
120061 120081 ). for which the barvonic contribution can be 
largely neglected and the mass density profile as derived 
from the RC "directly" traces that of the dark matter 
halo. The disagreement between an observed, apparently 
finite, central density and the expectation of a steeper 
density profile from the ACDM cosmological simulations 
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re mains an unsolved prob lem in galactic dynamics (but 
see iGovernato et aI]|201ClD . In fact, this cusp-core con- 
troversy has often been invoked as a major weakness of 
the current standard cosmological model of hierarchical 
growth of structures starting from a nearly homogeneous 
Universe, seeded with density fluctuations arising from a 
ACDM power spectrum. However, the impressive agree- 
ment between ACDM predictions and observations of the 
angular fluctuation spectrum of the Cosmic Microwave 
Background suggests that other explanations for the con- 
troversy must be sought, either in the dissipative and 
feedback physics of the baryons or in a more detailed 
analysis of halo density profiles from ACDM cosmologi- 
cal simulati ons. 

Recently, iNavarro et al.l ()2004D have proposed another 
model that fits the density profiles of halos in ACDM sim- 
ulations even better than th e NFW model. It was later 
re alized (iMerritt et al.ll2006f ) that the model advocated 
by INavarro et al.l (|2004[ ) had been previously introduced 
for the distribution of stel lar light and mass in galaxies 
(|Einastolll96"5l 119681 Il969l ). with density profile 



PE(r) = p-2 cxp ■ 



-2n 



r-2 



l/n 



- 1 



(3) 



where, again, r_2 is the radius where the density profile 
has a slope of —2 and p_2 is the local density at that ra- 
dius. While both the NFW and Iso models are described 
by two parameters, a characteristic scale and a charac- 
teristic density at that radius, the Einasto model involves 
a third parameter, n, the Einasto index, which describes 
the shape of the density profile. The Einasto mode l is the 
three -dimensional equivalent of the Sersic model (jSersid 
Il968| ) that provides an excellent fit to t he surface bright- 
ness profiles of elliptical galaxies (e.g. iKormendv et al.l 
120091) Fl 

The mass profile of th e Einasto halo 
(ICardone. Piedipalumbo fc Tortoral 120051 : 
iMamon fc Loka£ll2005D is written as 



r-2 



(4) 



where 7(3n, x) = e~^t^^~^dt is the incomplete gamma 
function. 

As illustrated in Fig. [T] (top panel), the net effect of 
increasing the halo shape index n at fixed characteristic 
density and radius is to increase and steepen the den- 
sity profile p(r) in the central part of the halo. This then 
leads to increase the inner slope and amplitude of the RC 
(bottom panel of Fig.[T|). Equation [J] also implies that an 
Einasto halo has a finite mass, whose physical feature is 
advantageous with respect to the divergent NFW and Iso 
mass profiles. At fixed characteristic Einasto scale den- 
sity and radius, the mass profile converges more rapidly 
to its virial mass for small indices than for larger ones. 

Because the inner RCs are expected to probe deeply 
inside the halos down to scales of 1 /300th of the virial 
radius, where Navarro and co-workers have found small 
but significant departures between the density profiles of 



5 Both the NFW model IILokas fc MamonI [200lD and halos m 
ACDM simulations ijMerritt et al.ll2005l ) are well-fit by an n = 3 
Sersic model (Sersic 1968). 
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Fig. 1. — Density profile {top) and rotation curve (bottom) for 
Einasto models of different indices. The characteristic radius r_2 = 
20 kpc and density p_2 = 1.5 X 10~^ M<3 pc~^ are the same for 
all models. Only the Einasto index n changes as indicated by 
values and lines of various colour and style. The dotted black lines 
correspond to models with n = 0.005 (bottom curve) and n = 0.1 
(upper curve). 



ACDM halos and the NFW model, it has been suggested 
(e.g. [Stoehr 200^ that the shallower inner slope of the 
Einasto model should reconcile the observed RCs with 
ACDM models. However, no extensive mass decomposi- 
tions of galactic RCs using Einasto halos have yet been 
carried out. Only the RC of the Andromeda gala xy has 
been modeled with the Einasto formula (Chemi n et al.l 
l2009f l. though with no real improvement with respect to 
the usual core and NFW halos as cause d by the pecu- 
liar sh ape of the M31 RC. Note also that iGraham et al.l 
(|200a compare the central dark matter densities of a 
sample of low surface brightness galaxies (|de Blokll2004[) 
to a family of Einasto halos and conclude that these 
match the data reasonably well. 

The present article aims to provide the first mass de- 
compositions of RCs with the Einasto model on an im- 
portant sample of spiral galaxies. The questions we wish 
to address are the following. Is the Einasto halo a good 
fit to the RCs of galaxies? Is it a better description of 
observations than the usual NFW and Iso halos? What 
does a typical galactic Einasto halo look like? 

We use a subsample from The HI Nearby Galaxy 
Survey (THINGS), as described in lde Blok etall (120081 
hereafter DOS) and in ^ and compare the quality of 
Einasto halo fits with those obtained using the NFW 
halo and pseudo- isothermal sphere (0. We investigate 
the statistical significance of possible improvements of 
the Einasto model O), determine the parameter space 
of Einasto halos for the THINGS sample (iJ5|) and com- 
pare the properties of galaxy-sized halos generated in 
numerical simulations to those derived from our sample 
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(Ijni)- We also describe a range of Einasto indices (con- 
strained two-parameter models) that provide better fits 
to RCs than the Iso and NEW models (SjT]). This ar- 
ticle is intended to be the first in a series that aims at 
investigating how well the Einasto formalism works in 
galactic dynamics and at quantifying the degree of cus- 
piness of dark matter hal os as a function of g alaxy mass. 
The WMAP3 cosmology (jSpergel et al.ll2007l ) is adopted 
throughout this work. 

2. THE THINGS SAMPLE 

The Hi observations of galaxies used in this study 
are from The Hi Nearby Galaxy Survey (THINGS, 
iWalter et all 120081 ). This survey consists of high- 
resolution 21-cm observations of thirty-four nearby 
(closer than 15 Mpc) spiral and irregular galaxies (Sb 
to Im). The galaxies were observed with the Very Large 
Array (VLA) in its B, C, and D configurations. The 
Hi data and the complete details of their reduction have 
been present ed in previous articles from the THINGS se- 
ries (see e.g. IWalter et al.ll2008f ). We analyze here the 
sample of seventeen rotationally dominated and undis- 
turbed galaxies presen ted in D08, which i s a subset of 
the whole sample from IWalter et all ()2008[) . 

3. MASS MODELING 

3.1. Generalities 

The reader is referred to D08 for the details of the 
kinematical analysis of the Hi velocity fields by tilted 
ring models to extract the rotation curves. D08 has also 
presented a complete dynamical analysis of those RCs 
with the NEW and Iso models. 

We have modeled the mass distribution using differ- 
ent RCs from those published in D08, as explained in 
N3.2I We have performed Levenberg-Marquardt non- 
linear least-squares fits to the RCs, taking into account 
the contributions of the g spherical stellar bulge, a 
stellar disk, and a spherical dark matter component. 

These contributions were taken from D08 for the stel- 
lar and gaseous RCs and from equation (U) for the 
halo. Basically the (atomic only) gaseous disk RCs come 
from THINGS Hi surface densities and the stellar RCs 
from 3.6 fini surface brightness profiles, available for all 
our galaxies from complementary obser vations with the 
Spitz er Nearby Galaxy Survey (SINGS. iKennicutt et al.l 
I2003f) . The advantage of using these near-infrared lu- 
minosities, is that they are little affected by internal 
or Galactic extinction. We have not considered mod- 
els with free mass-to-light ratios in combination with 
Einasto halos. Such an analysis will be presented in a 
future paper (Chemin et al., in preparation). Instead 
we have adopted the fixed mass-to-light ratios Ta g/^m 
of DOS, who derived T3.6^,„ from the 2MAS S J - 
color s, adopting T3 = 0.92 Tk^ — 0.05 (|0h et al.l 
12001 and logT^ = 1.43 ( J - K,) - 1.38, extrapol ated 
from IBell fc de Jond ([200lh . IBell fc de .TomI ([2001 de- 
rived stellar mass-to-light ratios assuming stellar popu- 
lations with a bursty star formation history with a diet- 
Salpetcr initi al mass f unctio n (IME). As in DOS we also 
consider the iKroupal ()200lD IME, which decreases the 
mass-to-light ratios (hence stellar masses) by 0.15 dex. A 
sech^(z/zo) law ()van der Kruit fc Searlel[l98l was used 
to derive the vertical distribution of the stellar disk, un- 
der the hypothesis that the disk scale-height zq is one 



fifth t he radial disk scale - length ()van der Kruit fc Searld 
fT98llKregel et al.l[200l) . 

None of our mass models take into account the contri- 
bution of molecular gas because the total gas surface den- 
sities ar e atomic gas dom inated for the majority of our 
sample (jLerov et al.ll2008l ) . However we briefiy report on 
the effects of molecular gas on Einasto halo parameters 
for the extreme case of NGC 6946 in gSH 

Following DOS the RCs of the galaxies NGC 2403 and 
NGC 3198 have been decomposed using either one sin- 
gle stellar disk component (labeled "NGC2403d" and 
"NGC319Sd") or two stellar components. Also, for 
NGC7793 we present the analysis of the whole RC, as 
well as that of its rising part only ("NGC7793s" ). For 
each assumed IMF, fits have been performed using fixed 
and free Einasto indices for sake of comparison with 
cosmological simulations, as will be presented in Sec- 
tion [6l Furt her technical details of the fitting procedure 
are given m IChemin et all ([20091) . 

3.2. The influence of RC sampling on mass models 

The RCs we fit here are sampled differently from those 
analyzed in DOS, where the Hi rotation velocities used for 
the dynamical analysis are not totally independent. DOS 
used a sampling of two data points per spatial resolution 
element. As a consequence they did not directly compare 
the values obtained for the core and cusp models but 
instead studied the differences in x^- 

Since a direct comparison of the reduced values for 
the different two- and three-parameter models is essen- 
tial to identify the "best" halo model, one has first to in- 
vestigate whether any correlation between spatially adja- 
cent rotation velocities could alter the results of the mass 
modeling. For this purpose, we have deri ved mass mod- 
els of several galactic RCs from DOS and IChemin et al.l 
(|200a using three different samplings, assuming Einasto, 
NEW and Iso halos. 

A first modeling exercise (SI) was done with RCs sam- 
pled with two points per beam (as i n DOS), a second one 
(S2) with one point per beam (as in IChemin et al.ll2009| ) 
and a third one (S3) with one point every two beams. 
The last two samplings provide independent data points, 
as opposed to the fir st one. The velocity u ncertainties are 
defined as in DOS or I Chemin et al.l ()2009D . taking into ac- 
count the (generally dominant) kinematical asymmetry 
observed between the receding the approaching halves of 
each galactic disk. 

We observe the following results (which are strictly 
speaking only valid for the range of samplings explored 
here and should not be further generalized) : (i) the pa- 
rameter space of the dark matter halo models does not 
vary with the spatial sampling of their RC, (ii) the rel- 
ative fit quality of the three halo models does not vary 
with resolution, and (iii) the x^ values increase from SI 
to S2 and then decrease for S3, which generally exhibits 
the smallest x^ ■ To illustrate, the fits of the NEW, Iso 
and Einasto halos to the RC of the Andromeda galaxy 
lead to a reduced ^ for SI (sampling of ~200 pc), 
21 for S2 ( a sam pling of ~400 pc, as in Tab. 5 
of IChemin et all ()2009D . for their "hybrid" model) and 
~ S for S3 (samphng of ~S00 pc). For this galaxy, the 
fit parameters of the NEW halo, the characteristic veloc- 
ity ^200 and the halo concentration c, are V200 = (145±2) 
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Fig. 2. — Comparison between reduced from mass models fit 
to THINGS RCs from DOS and from this study, with the Iso (Open 
diamonds) and NFW {closed diamonds), using the Kroupa (green) 
and dict-Salpeter (purple) IMFs. 

km and c = 21 ± 1 for SI, ^200 = (146 ± 4) km s'^ 
and c = 20 ± 2 for S2, ^200 = (147 ± 5) km s^^ and 
c = 21 ± 2 for S3, indeed independent of the sampling. 

The non-independent rotation velocities in SI are thus 
responsible for the smaller compared to S2, whereas 
the larger sampling interval explains the smaller for 
model S3. The risk of using large intervals (like S3) is, of 
course, that one starts loosing small-scale details in the 
rotation curve and decreases the velocity gradient in the 
inner parts of the galaxies. 

We therefore conclude that the optimum sampling of 
RCs for mass models of galaxies is one which avoids 
non-independent velocities generating artificially small 
X^, and which preserves an appropriate spatial resolu- 
tion that does not smooth out small-scale details in the 
RC, especially in the inner regions, which are the most 
sensitive to the differences between the Iso, NFW and 
Einasto models. We have thus modeled the mass distri- 
butions of the THINGS galaxies from RCs with the S2 
sampling, i.e., with galactocentric radii separated by one 
synthesized beam size. 

Note that the conclusions given in DOS still hold: the 
halo parameters we have fit are similar to theirs, the 
Iso model usually provides better results than the NFW 
model. Only the values of the reduced have increased 
when going from non-independent to independent veloc- 
ities, as shown in Figure [5J Finally, because the analysis 
of individual differences between the NFW and Iso halos 
is beyond the scope of this article, we refer the reader 
to Tab. [1] and [2] for the resulting halo parameters and to 
DOS for further discussion. 

3.3. The Einasto model fits 

Figures[3]and|3]show the decomposition of the RCs into 
baryonic and dark matter components for the Einasto 
model. The fit parameters and their associated la error, 
as well as the reduced chi-squared values, x^, are given 
in Tabs [T] and [2] for the Einasto, NFW and Iso halos. 
Figure [5] displays the values of Einasto Xr f^nd compares 
them with NFW and Iso Xr values. 

The most important results from the analysis of the 
fits are: 



• The quality of the fits is very good for the major- 
ity of the galaxies. S0% (16 out of 20) of the fits 
done with a Kroupa IMF have Xr < 1-5- This frac- 
tion falls to 60% (12 out of 20) when using a diet- 
Salpcter IMF. The Einasto model is thus highly 
constrained (low x^) owing to the large numbers of 
degrees of freedom of each fit. 

• Most of the fits (70%) with a Kroupa IMF are bet- 
ter (i.e., have lower Xr) than those with a diet- 
Salpeter IMF. That trend is in good agreement 
with the results obtained for the NFW and Iso 
halos (see DOS and our current study in Tabs. [T] 
and m . Both our analysis and that of DOS show 
how so me mass models done under the prescrip- 
tions of iBell fc de Jongj ()200ir ) with the assump- 
tion of a diet-Salpeter IMF can be unphysical due 
to modelled velocities exceeding observed ones in 
the inner disk regions (for, e.g., NGC 3521, NGC 
5055, NGC 6946 or NGC 7331). It explains why 
the Einasto halo properties deduced from the diet- 
Salpter IMF often present scattered distributions 
and/or large errors, as will be shown in following 
sections. It is thus tempting to reject that IMF hy- 
pothesis. It is the reason why we focus more on the 
results obtained with the Kroupa IMF hypothesis 
hereafter, except where mentioned. 

• The Einasto models give better results than the 
Iso and NFW models, irrespective of the IMF: 60% 
(S0%) of the Einasto halos fit better than Iso ha- 
los for the dict-Salpeter IMF (Kroupa IMF, respec- 
tively); all of them fit better than the NFW cusp 
model, irrespective of the IMF. This result is ex- 
pected owing to the effect of the third modeling 
parameter (the Einasto index) in addition to the 
usual characteristic scale halo density and size. We 
analyse the statistical significance of that modeling 
improvement in |21 Furthermore it is likely that the 
exponential-like decrease of the dark matter volu- 
mic density described by Eq. |3] contributes to the 
fit improvement as well, as will be shown in fJZl 

4. THE EINASTO MODEL AS THE PREFERRED MODEL 

In this section, we investigate the significance of the 
improvement of the RC fits with the Einasto halo over 
those with the Iso and NFW models. 

We first analyze the number of fits that are statistically 
significant with the Einasto model and that were not 
significant with any of the two-parameter Iso and NFW 
models. By "statistically significant" we mean that a 
model and the observations follow the same distribution, 
i.e., that there is less than 5% probability of obtaining as 
good a fit by chance according to chi-squared statistics. 
For the Kroupa IMF, the fraction of the fits that were not 
significant with the Iso model, but are with the Einasto 
halo consists of 30% of the sample (6 out of 20 fits). This 
fraction becomes 25% when going from NFW to Einasto. 
For the diet-Salpeter IMF, only two galaxies (10%) have 
gone from a not significant result with the Iso model to 
significant with the Einasto halo, while it is 20% from 
NFW to Einasto (4 galaxies). 

Though these numbers do not represent the majority 
of the sample, they arc not negligible. As a comparison. 
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TABLE 1 

Fit parameters of the Einasto halo to THINGS galaxies for a free index and fixed mass-to-light ratios derived using a 

diet-Salpeter IMF. 
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4.1 it 0.1 


0.4 


27.6 it 2.4 


1.3 it 0.1 


NGC5055 


7.7 


0.7 it 0.2 


32.2 it 4.3 


0.1 it 0.2 


10.7 


234.6 it 26.1 


0.9 it 0.1 


8.1 


0.9 it 0.3 


44.7 it 30.2 


NGC6946 


1.3 


4.6 it 2.9 


14.9 it 8.3 


2x10-3 ± 0.3 


2.8 


307.4 it 37.5 


1.9 it 0.3 


1.5 


5.3 it 0.5 


21.4 it 5.8 


NGC7331 


3.1 


1.0 it 4.8x10^ 


> 1000 


0.2 it 3x10^5 


4.3 


275.0 it 131.7 


1.1 it 0.9 


3.0 


1.6 it 0.4 


113.3 it 519.4 


NGC7793 


1.7 


22.7 it 1.9 


3.7 it 0.1 


0.3 it 0.1 


4.4 


141.4 it 8.1 


6.6 it 0.2 


3.1 


78.2 it 11.2 


1.9 it 0.2 


NGC7793S 


1.5 


< 10-3 


> 1000 


9.3 it 42.3 


4.8 


181.0 it 20.8 


5.6 it 0.4 


2.2 


51.8 it 6.8 


3.5 it 0.7 



Note. — Columns (l)-(5)-(8) Reduced . Columns (2)-(9) Characteristic density. Columns (3)-(10) Characteristic radius. Column (4) Einasto 
index. Columns (6)-(7) NEW halo characteristic velocity and concentration. 



TABLE 2 

Same as in Tab.[1]but with mass-to-light ratios derived using a Kroupa IMF. 
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NEW 
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Galaxy 
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n 
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r 


V200 
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X 


2 

r 


PO 


rc 










(10-3 Mq pc-3) 


(kpc) 








(km s— -'^) 








(10-3 Mq pc-3) 


(kpc) 


NGC925 





6 


4.5 it 0.5 


8.9 it 0.6 


0.1 it 0.1 


3 


4 


171.4 it 15.2 


2.2 it 0.1 


1 


1 


6.0 it 0.8 


9.6 


it 1.8 


NGC2366 





2 


5.7 it 1.1 


3.2 it 0.3 


1.1 it 0.2 


1 


2 


75.4 it 16.3 


4.2 it 0.8 





2 


40.1 it 4.2 


1.3 


it 0.1 


NGC2403d 





6 


1.5 it 0.4 


15.6 it 2.0 


5.6 it 0.5 





6 


101.2 it 1.0 


12.6 it 0.3 


1 


1 


154.0 it 10.8 


1.5 


it 0.1 


NGC2403 





6 


1.7 it 0.4 


14.6 it 1.8 


5.3 it 0.5 





6 


102.0 it 0.9 


12.4 it 0.2 


1 





147.0 it 9.9 


1.5 


it 0.1 


NGC2841 





2 


10.7 it 1.2 


11.2 it 0.6 


12.4 it 1.5 





6 


172.0 it 1.4 


25.2 it 0.6 





2 


> 1000 


0.6 


it 0.1 


NGC2903 





3 


72.4 it 19.4 


3.2 it 0.4 


9.6 it 2.1 





4 


111.2 it 0.8 


36.2 it 1.0 


1 


2 


> 1000 


0.0 


it 0.1 


NGC2976 





5 


0.0 it 0.6 


> 1000 


4.9 it 77.8 


2 


3 


108.6 it 5.9 


4.5 it 0.2 





5 


36.6 it 4.7 


4.5 


it 2.8 


NGC3031 


3 


1 


30.0 it 3.0 


4.7 it 0.2 


0.7 it 0.2 


3 


6 


94.3 it 5.5 


26.8 it 3.6 


3 


9 


740.0 it 434.0 


0.8 


it 0.3 


NGC3198d 





8 


2.2 it 0.2 


13.1 it 0.6 


2.1 it 0.3 


1 


4 


109.3 it 2.5 


8.9 it 0.6 





8 


47.1 it 5.9 


2.7 


it 0.2 


NGC3198 


1 


5 


2.2 it 0.3 


13.0 it 0.7 


2.0 it 0.3 


2 


2 


109.9 it 3.2 


8.7 it 0.7 


1 


5 


44.3 it 7.4 


2.8 


it 0.3 


IC2574 





2 


0.7 it 0.3 


15.1 it 5.0 


1.1 it 0.3 


3 


4 


94.6 it 4.3 


2.2 it 0.2 





2 


4.9 it 0.2 


6.3 


it 0.3 


NGC3521 


4 


9 


6.1 it 2.5 


9.4 it 1.6 


1.2 it 0.9 


5 


8 


121.7 it 22.1 


9.7 it 3.1 


4 


8 


75.6 it 45.8 


2.5 


it 0.9 


NGC3621 





6 


0.3 it 0.2 


39.9 it 12.3 


6.4 it 1.0 





6 


119.8 it 1.9 


7.9 it 0.3 


1 


3 


49.3 it 4.2 


2.8 


it 0.1 


NGC4736 


1 


5 


50.8 it 15.5 


1.8 it 0.3 


1.9 it 0.8 


1 


5 


44.1 it 2.3 


43.8 it 9.0 


1 


8 


> 1000 


0.1 


it 0.1 


DD0154 





3 


1.5 it 0.3 


5.9 it 0.8 


2.2 it 0.3 


1 


4 


91.3 it 4.0 


2.7 it 0.1 





4 


28.5 it 2.5 


1.3 


it 0.1 


NGC5055 





8 


1.3 it 0.1 


22.6 it 0.6 


0.5 it 0.1 


1 


6 


234.4 it 22.5 


1.9 it 0.3 


1 





4.8 it 0.5 


11.7 it 0.9 


NGC6946 


1 





1.8 it 1.2 


18.7 it 7.4 


2.6 it 0.9 


1 





187.9 it 17.0 


6.1 it 0.7 


1 





44.2 it 4.1 


3.7 it 0.2 


NGC7331 





3 


0.0 it 0.0 


> 1000 


13.1 it 19.7 





3 


195.6 it 15.3 


5.2 it 0.6 





4 


26.1 it 3.6 


5.2 it 0.5 


NGC7793 


2 


8 


19.8 it 3.7 


3.6 it 0.3 


0.9 it 0.2 


3 


9 


114.4 it 10.2 


9.2 it 0.8 


3 


7 


128.0 it 17.6 


1.5 it 0.1 


NGC7793S 


1 


7 


< 10-3 


> 1000 


12.7 it 32.9 


2 


9 


191.9 it 9.7 


6.1 it 0.4 


4 


1 


96.7 it 15.1 


2.0 it 0.3 



15% (20%) of the sample acquires significant fits when 
passing from NFW to Iso for the Kroupa (respectively 
diet-Salpeter) IMF. Therefore, the improvement of fit- 
ting the Hi RCs by the Einasto model instead of the Iso 
or NFW models is equal to or better than the improve- 
ment implied when switching from NFW to Iso. Con- 
versely, the Einasto model does not diminish the signifi- 
cance of the fits : no fit has turned from being significant 
with Iso/NFW to not significant with the Einasto halo, 
regardless of the IMF. 

However, it is not surprising that the three-parameter 
Einasto model is better able to reproduce the RCs than 



the two-parameter Iso and NFW models. Because nei- 
ther Iso nor NFW models are nestecj^ into the Einasto 
family of models, we cannot simply perform an _F-test 
to check for the significance of the improvement of the 
fits with the three-parameter Einasto model. Instead, we 
verify the significance of this improvement by perform- 
ing Akaike i nformation criterion tests (hereafter, AIC, 
IAkaikelfT97l . which are designed to compare non-nested 

" A model with k parameters is considered nested in a more 
complex model with e.g., k + I parameters when fixing one of the 
k + 1 parameters allows one to reproduce the simplest k parameter 
model. 
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Fig. 3. — Mass models for the Einasto halo derived from the Hi RC of THINGS galaxies. Displayed results arc for fixed stellar mass- 
to-light ratios derived from stellar population synthesis with a diet-Salpeter initial mass function and whose values are given in DOS. A 
red dashcd-dottcd (dotted) line is for the stellar disk (bulge, respectively) components, a blue dashed line for the atomic gas, a green 
dashcd-dottcd line for the dark halo and an orange thick line for the overall model. 



models of different numbers of parameters. The goal of 
an AIC test is not to statistically reject one model with 
respect to another, but rather to find the model that 
is more likely to be correct. The criterion is expressed 
as AIC = + 2-/V, where N is the number of parame- 
ters of the model. An Einasto model is considered more 
likely to be correct than the Iso or NFW model when 
AICEinasto < AICi,o/NFW- Wc find that 65% (60%) of 
the fits with the Einasto model are more likely to be cor- 
rect than those with the ISO model assuming a Kroupa 
(dict-Salpeter) IMF. These percentages increase to 95% 
when comparing with the NFW model, independent of 
the assumed IMF. 

In summary, the Einasto model provides more fits that 
are statistically significant and is more likely to be cor- 



rect than the two other models. For these reasons, we 
conclude that fitting RCs of the current THINGS sample 
with an Einasto halo is a significant improvement with 
respect to the NFW and Iso models. 

5. DYNAMICAL PROPERTIES OF EINASTO HALOS 

5.1. Non-universality of the halo 

Figures [6| and [7| show the fit parameters for the diet- 
Salpeter and Kroupa IMFs, respectively. Figure [8] shows 
the distribution of Einasto indices. Based on these dis- 
tributions, we can distinguish several families of halos. 
A first family of halos has well constrained parameters 
with large indices {n > 4). We refer to this family as 
cuspy Einasto halos hereafter. Indeed it corresponds to 
halos whose indices are comparable with the index of 
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Fig. 4. — Same as Fig. [3] but with a Kroupa initial mass function. 
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Fig. 5. — Comparison of reduced for the Einasto halo with the 
NFW and Iso halos (respectively filled and open symbols), using 
the diet-Salpoter [left) and Kroupa (right) IMFs. 



a typical galaxy-sized halo from numerical simulations 



(5 < n, 



.simiilatRH 



< 7, iNavarro et all l2004t iMerritt et al.l 



[2506t IGraham et alll2006l : INavarro et all [20101) . Their 
innermost density slopes are steeper than —0.8, hence 
consistent with cos mological slopes (a = —0.9 ± 0.1, 
INavarro et al.l I2010D . A second family of galaxies has 
well-constrained parameters with low index values (0.1 < 
n < 4). It contains the majority of the galaxies. A 
third halo family has very low indices (n < 0.1). Most 
of those fits occur with the diet-Salpeter IMF. Finally 
a fourth family of halos has unrealistic, extremely large 
scale radii, with large associated uncertainties (the ris- 
ing part of the RC of NGC 7793, NGC 7331, NGC 2976). 
Those fits are extremely degenerate, which is very likely 
caused by an almost complete dominance of the baryonic 
material over the dark component, at least for NGC 2976 
and NGC 7331. 
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Fig. 6. — Parameter space [p_2i ''-2i of Einasto halos for the THINGS sample (filled diamonds). Fixed stellar mass-to-light ratios 
arc derived from stellar population synthesis using a diet-Salpctcr initial mass function. Colored sym bols show the parameters of dwarf- 
(upward r ed triangles) and gal axy- (downward green triangles) halos modeled in ACDM simulations ijNavarro et al. 20(53). Galaxy-sized 
halos from lTissera et al.l 1120101 ) are displayed with orange circles. A dashed line is a power law fit to the relationship linking p_2 and r—2- 
log(p-2) = (—1.61 ± 0.07) log(r_2) — (0.81 ± 0.06). The degenerate fits with log r-2 > 3 and log p-2 > —5 arc not drawn for clarity. 



Figures|6]and[7]show a clear relation linking the charac- 
teristic density to the radius. Small halos are denser than 
large ones. A power-law fit to the observed relationship 
gives p_2 oc rZl'^^"'^ for both IMFs. This implies that 
the surface density of dark matter I]_2 = p_2?'-2 derived 

at r-2 scales as S_2 oc rZj'^ oc Those relations 

have been determined without taking account the fourth 
halo family fits and by averaging the fit parameters of 
galaxies that have duplicated RCs but a different stellar 
decomposition (NGC 2403 and NGC 3198). It is worth- 
while to note that the correlation we find between the 
scale density and radius of the Einasto model differs from 
the relations found between the scale density and radius 
of two-parameter NFW and Iso m odels. For h i stance 
iKormendv fc FreemanI ()2004[) and iSpano et all (I2008D 
have derived a relation po oc 1/ro (but see iBarnes et alJ 
120041 ) ■ where the scale density po and and t-q are the 
characteristic density and radius of two-parameter NFW 
and Iso models. As a consequence the dark matter sur- 
face density Eg = po r emains roughly co nstant at 
the core radius of the halo (IKormendv fc Freem an 200j; 
ISpano et al1l2008HDonato et al.ll2009[) . This trend is not 
observed with the Einasto halo. In this context, it is in- 
teresting to note that we have also fit the two-parameter 
Iso model to the current sample (Iso column in Tabs. |T] 
an d 121 ■ The results do not imp l y a similar r e lation as 
in IKormendv fc Freema"3 (|2004l ) . (Spano et all (I2008D or 
IDonato et al.l (j2009D . but agree more with IBarnes et al.l 
(|2004D . We indeed find inner surface densities for the 
core halo that depend on the scale parameters. A lin- 
ear fit yields po oc r^^-^^° °^ for the diet-Salpeter IMF 
and pq oc 7--1-5±o.04 Kroup a IMF. It may seem 

odd that these results disagree with lDonato et alT()2009l ) 
because part of their sample consists of the current 
THINGS galaxies (though using the DOS RCs). A pos- 
sible explanation of the differ ence could be the choice 
of stellar mass-to-light ratios (jPonato et al.l use best fit 
mass-to-light ratios) , the use of other kinematical tracers 
(mostly the ionized gas) as well as the choice of the outer 
slope of the den sity profile of dark matter {p oc in 
ISpano et al.l and IDonato et al.l . while our Iso model im- 
plies p oc r~^). A detailed investigation of both analy- 



sis methods is beyond the scope of this article. Clearly 
though, the implied (non-)constancy of the dark matter 
surface density warrants further investigation. 

As for the Einasto index, we note the more scat- 
tered distribution and the presence of more halos hav- 
ing a small index for the dict-Salpeter IMF than for the 
Kroupa IMF. This is caused by a larger contribution of 
stellar baryons to the RCs. The effect of a (very) low 
index (illustrated in Fig. [T] as declining velocities in the 
outer parts of the dark matter halo) are clearly visible for 
NGC 4736, NGC 6946 (Fig. El diet-Salpetcr IMF only) 
and NGC 925 (Figs. O and IH). 

No obvious correlation is found between the Einasto 
index and scale density or radius. Using halos with 
n < 4 (second and third halo families) we derive an 
average index of n = 0.8 ± 0.3 with the diet-Salpeter 
IMF and fi = 1.3 ± 0.2 with the Kroupa IMF (us- 
ing averaged indices for NGC 2403/NGC 2403d and 
NGC 3198/NGC 3198d). Note that the amplitude and 
scatter become more important if the first halo family is 
also taken into account, with n = 2.1±0.8 (diet-Salpeter) 
and n = 3.2 ± 0.9 (Kroupa). 

The cuspiest Einasto halos are typically found in galax- 
ies where the dark matter dominates the visible matter 
at almost all radii (as seen in Figs. [3] and |4]), coupled with 
a relatively extended rotation curve and mass distribu- 
tion {r > 15 kpc) with rotation velocities larger than 120 
km s^^ at large galactocentric radius. Note also they 
correspond to galaxies whose RCs start at r > 3 kpc. 
The missing information at r < 3 kpc might affect our 
conclusions for those galaxies. 

Perhaps the most important result from the analysis of 
the indices is that no single value or well-defined relation 
can reproduce all types of halos. The index is also a 
complex function of the halo mass (see H5.3^ . We can 
thus conclude that no obvious "universal" Einasto index, 
hence Einasto halo, can be deduced from the current 
galaxy sample. 

5.2. Intermediate class halo density and density slope 

profiles 

Figure |9] displays the radial profiles of the mass vol- 
ume density and the radial profiles of the logarithmic 
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TABLE 3 

Logarithmic slope a of the dark matter density 

PROFILE DERIVED AT log(r7r_2) = —1.5. 



-3 -2-10 1 
log(Einasto index) 

Fig. 8. — Distribution of Einasto indices. Results for the diet- 
Salpeter and Kroupa IMFs are shown by blue and red histograms, 
respectively. The range of Einasto indices from ACDM simulations 
with and without the physics of baryons (5 < n < 10) is displayed 
by dashed lines. 

density slope of the THINGS sample. For clarity, only 
results for the Kroupa IMF are represented, excluding 
the fourth halo family. Our conclusions are similar for 
the diet-Salpeter IMF. Also shown are the radial profiles 
of generic galaxy-sized Iso and NFW halos with compara- 
ble densities, giving a halo RC with a velocity amplitude 
of ~ 200 km s~^ at a typical radius of 25 kpc. Other 
dashed colored curves correspond to Einasto halos with 
parameters derived from cosmological simulations (see ^ 
for details). 

Figure [5] (top panel) shows that the density profiles 
typically decrease very smoothly in the inner regions. 
Further out, the profiles become very steep and the den- 
sities very low past the scale radius r_2. It can also be 
seen that most of the halos derived from the observations 
are less dense in their inner parts than any of the simu- 
lated galaxy-sized halos. All of this is explained by the 
fact that THINGS galaxies generally have small Einasto 
indices whereas simulated halos have large index values 
due to their cuspy nature (see ^ for a complete discus- 
sion) . 

The bottom panel of Fig.[S]shows that on the one hand, 
the logarithmic slopes in the innermost regions are more 
reminiscent of the slope of a Iso halo than that of cos- 
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NGC7793 
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mological Einasto or NFW cusps, for the majority of 
the sample. On the other hand, because of the progres- 
sively steepening nature implied by the Einasto model 
(Eq. [3]) , the observed slope profiles become similar to a 
NFW cusp at larger radii. The slopes are even larger 
than those of any cosmological halos beyond r_2, as a 
result of the small amplitude of the observed Einasto 
indices. A density slope profile of a typical observed 
Einasto halo (log(p_2 -A^opc"'^) ^ —2, r_2 ~ 10 kpc, 
n 1.3, as fit with the Kroupa IMF) thus represents an 
intermediate case between a pseudo-isothermal core halo 
and cosmological cusps inside r_2. As a consequence of 
these results, we estimate that half of the total mass of 
such halos is contained within about 15 kpc, and the 
mass profile has rapidly converged towards, e.g. 95% its 
total mass at r ~ 40 kpc only (or four times the average 
scale radius r_2). 

We report the inner slopes a = dlog(p)/dlog(r) of the 
dark matter density profiles in Tab. [3] For sake of uni- 
formity with cosmological simulations whose inner slopes 
are derived at a fraction of the Einasto radius r_2 (ba- 
sically —2 < log(r/r_2) < —1, Navarro et al. 2010), we 
have derived the inner slopes of the THINGS sample at 
log(r/r_2) = —1.5 (or r = 0.03r_2), discarding all fits 
of the fourth halo family. The mean logarithmic slope of 
the density profiles is a = — 0.3±0.1 for the diet-Salpeter 
IMF and a = -0.5 ± 0.1 for the Kroupa IMF. Those 
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numbers become a = —0.1 ±0.1 and a — —0.2 ± 0.1 (re- 
spectively) when the cuspiest halos of the sample {n > 4) 
are removed from the distribution. Note that the mean 
inner density slope is — 1.3 ± 0.2 for those cuspiest halos 
only. For comparison, the innermost slope in the Aquar- 
ius pure-ACDM simulation (at r_2/100) is -0.9 ± 0.1, 
which is shallower than the cuspiest of our halos, but 
steeper than most of them. 




Fig. 9. — Dark matter density profiles (top) and density slope 
profiles (bottom) for the fit Einasto halos (black, solid lines). Only 
results for the Kroupa IMF are shown for clarity. For each halo 
the radial range of the observations is highlighted by a thick line. 
Colored dashed-dotted lines represent Einasto halos from ACDM 
simulations whose parameters are shown in Fig.|6] A purple dashed 
line represents a generic galaxy-sized NFW halo, a magenta dotted 
line a generic galaxy-sized Iso halo. 



As a reminder, an almost constant dark matter den- 
sity has been found for a la rge sample of dark mat - 
tcr dominated galaxies (e.g. Ide Blok fc Bosmal I2002D . 
These authors have found a mean inner density slope 
a ^ —0.2, as measured "directly" from optical spec- 
troscopy by inverting the RCs into volume densities 
and extracting the slope at the innermost data point 
of the observations. We do not know yet the aver- 
age slope that would be fit with the Einasto profile 
for similar low surface brightness objects. The current 
THINGS sample consists mostly of g alaxie s more mas- 
sive than those in Ide Blok fc Bosmal (120021 ) (stellar and 
total masses), despite the presence of a few low surface 
density objects (among which NGC 2366 is common with 
Ide Blok fc Bosmall2002f ). Note, however, that both sam- 
ples show comparable average inner density slopes of 
dark matter. At first it may seem surprising because 



the quantity of baryons strongly differs in those samples. 
That particularity will be investigated in a future paper 
of our series. 

5.3. Halo concentration versus Einasto index and mass 

We have derived the virial radius r2oo and halo concen- 
tration C200 = ?'20o/''-2 as the radius of a sphere of mean 
density /02oo which equals two hundred times the criti- 
cal density for closure of the Universe, pcrit = 3i?Q/87rG'. 
Wc derive the (halo) mass M200 at r2oo using Eq.g) Fig- 
ure [TO] compares the virial masses, halo concentrations 
and indices. Concentrations and virial masses are re- 
ported in Tab. |4l For clarity, the fourth halo family has 
been omitted in Fig. [TO] and Tab. H] 

Adopting the Kroupa stellar mass scaling, it is ob- 
served that the Einasto index increases with the halo 
mass for halos more massive than M200 ^ 2 x 10^^ 
A4q while its distribution is more scattered for less mas- 
sive halos, showing an almost flat part. That curved 
trend is less evident with the diet-Salpeter IMF results. 
Though the mass range spanned by the observations is 
relatively small in comparison to what can be probed by 
numerical models, the complex relation between Einasto 
index and mass illustrates that no "universal" galactic 
dark matter halo can be identified from these observa- 
tions: the Einasto index cannot be simply deduced from 
the halo mass, nor does a simple scaling law allows one 
to scale the Einasto index/mass profile to those of a halo 
with a different mass. Our observational result thus con- 
firms the non-unive rsality of Einasto ha los seen in nu- 
merical simulations ([Navarro et al.|[2"010[) . 

The concentration is not correlated with the Einasto 
index (right-hand panel), but tends to decrease with the 
mass (middle panel). The scatter of the trend is how- 
ever significant. Halo fits with the diet-Salpeter IMF are 
less concentrated than halo fits with the Kroupa IMF, as 
caused by the more important contribution of the stellar 
component. As a comparison, the expected dependen- 
cies of the Einasto and NFW concentrations on the halo 
virial mas s from c osmological simulations ([Neto et al.l 
[2007: M accio et all [2008: Gao et al. 200^ 1 are displayed 
as well. These authors have found concentrations basi- 
cally scaling with the total mass as C200 oc Mj^g^. If one 
omits the two outliers which are the most concentr ated 
and cuspie s t halo s, the relations of lNeto et al.l (I2007D and 
IGao et ahl (|2008f ) are in good agreement with the seven 
halos at the upper end of the C200 ~ .^200 trend for the 
Kroupa IMF, and thus above most of concentrations de- 
rived with the diet- Salpeter IMF. The NFW relation of 
iMaccio et all (|2008l ) is more consistent with the lower end 
of the remaining halos at the lower end of the concentra- 
tions. Notice however that the comparison is of limited 
value as the halo masses can be measured directly from 
the numerical simulations whereas we can only provide 
estimates for the observations. 

We finally note that possible correlations between the 
properties of the Einasto halos with those of the baryonic 
matter (absolute magnitude, baryonic gas masses, disk 
and bulge characteristic scalelengths, etc.) will be fully 
studied into a forthcoming paper of our series (Chemin et 
al., in preparation). A preliminary analysis nevertheless 
shows Einasto indices scaling with (total) stellar masses, 
but barely with e.g. the bulge-to-disk luminosity ratio. 
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Fig. 10. — Virial masses M20O1 concentrations C200 and indices of Einasto halos. Open (filled) symbols are for the diet-Salpeter (Kroupa, 
respectively) IMF. In the middle panel a green dashed-dotted line re presents the dep endency of the Einasto halo concentrations on the 
mass log{c20o) = 2.646 — 0.138 log(Af2()o) from numerical simulations iCao et al.ll200a . concentratio ns derived with fixe d Einasto indices), 
a blue dotted the dependency of the NFW halo concentrations log(c2oo) = 1-765 — 0.08 3 logfAf^nn) l|Maccio et al.ir2008l ) and a red dashed 
lino the dependency of the NFW halo concentrations log(c20o) = 2.121 — 0.1 log(M2oo) flNeto et al.ll20071) . Halos with log(n) ~ —3 are not 
displayed for clarity reasons. 



5.4. The influence of total gas densities on Einasto 
parameters 

The contribution of the molecular gas has been ne- 
glected in our study. This is justified since the total gas 
surface densities are dominated by those of the atomic 
gas for most of the sample. Even when molecules dom- 
inate the total gas density in the innermost regions of a 
few galaxies from our s ample (NGC 3521 NGC 4736, 
NGC 5055, NGC 6946: ILerov ^TeII 12003) . the stellar 
surface densities completely dominate the mass budget 
of the baryons in those regions. They are generally 
more than ten times larger than the molecular densi- 
ties, with the exception of NGC 6946 where the ratio 
of the surface density of molecules to that of the stars 
sometimes reaches a factor of 1/3 in its central kilopar- 
sec (jLerov et al . 2008)13 We have thus fit Einasto mod- 
els by adding the molecular gas component to the mass 
distributi on for NGC 6946 o nly, with H2 surface densi- 
ties from iLerov et al.l ()2008f ). The top panel of Fig. [n] 
shows the results of that mass distribution model. As 
a result, the quality of fit remains unchanged with the 
Kroupa IMF {Xr ~ 1-0), as well as p-2 and r_2 within 
the uncertainties, with p_2 = (2.7 ± 1.0) x 10^"^ Mq 
pc~^ and r_2 = (15.9 ± 3.7) kpc. The Einasto index be- 
comes n — 1.2±0.4, which is about half the value derived 
with no molecular gas {n = 2.6 ± 0.9). Such a difference 
of halo cuspiness is clearly seen in the bottom panel of 
Fig. [TT] (green curve), where the difference of the Einasto 
halo RCs reach 20 km s^^ at 3 kpc. Though NGC 6946 
is probably an extreme case where the influence of the 
molecular gas is maximum on the halo shape, it is likely 
that Einasto indices would be even lower had we included 
the molecules in our fits. We also think the total gas 
component could play a more important contribution to 
the total baryonic dynamical budget than in our current 
study if models were done under the assumption of free 
stellar mass-to-light ratios. However such an hypothesis 
of Ts.e^m is beyond the scope of the article and we defer a 

^ Since the stel lar masses we use are larger than those given by 
ILerov et al.l II2008I ). as explained in D08, that ratio of molecular gas 
to stell ar surface densities is even smaller in our current analysis 
than in ILerov et"all pOOSl V 



complete analysis of the influence of the total/molecular 
gas contribution onto Einasto models with free mass-to- 
light ratios to our forthcoming article (Chemin et al., in 
preparation) . 

6. STRENGTHS AND WEAKNESSES OF EINASTO ACDM 

HALOS 

We now compare more precisely the parameter space 
of observed Einasto halos with those derived from 



NGC6946, n = 1 .2 (with molecular gas) 
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Fig. 11. — Mass model of NGC 6946 with the molecular gas con- 
tribution ( top panel). The molecular gas RC is shown by a black 
dashed line. All other stellar and gas curves are from Fig. |4] The 
bottom panel displays the differences between RCs of the total 
models of NGC 6946 done with and without the molecular gas (or- 
ange solid line, in absolute values), and between the corresponding 
Einasto halo RCs (green dashed-dotted line). 
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TABLE 4 

Einasto halo concentration C200 and virial mass M200 



Galaxy 



log(c20o) 
(1) 



Clog 

(2) 



log{M200) 

(3) 



log(c200) 
(4) 



(5) 



log(M200 

(6) 



pushing the index towards larger values (n — > 10) than 
for the pure dark matter case (n ^ 6), while the ob- 
servations go to an opposite way {n — J> 1 — 2). Note 
ffi ough recent hydrodynam ical cosmological simulations 
by IGovernato et al.l (j2010D with strong SN feedback cou- 



pled with a high threshold for star formation lead to dark 
ma tter halos without cusps. As theoretically explained 
bv iPontzen fc Governatol (|201lD . episodic SN events in 
dwarf galaxies alter the dark matter distribution in non- 
adiabatic ways leading to cored profiles. 

Only a few galaxies have a parameter space in rough 
agreement with simulations within the quoted uncertain- 
ties, though with some subtle differences. For instance, 
NGC 2903 can be associated with a dwarf-sized simu- 
lated halo in the index-radius graph, but with a halo 
more massive than that of a dwarf galaxy (log(p_2) ~ —1 
A4q pc~'^) in the index-density parameter space. For 
NGC 3621 any (dis)agreement depends critically on the 
Note. — Columns (l)-(2)-(3) arc for the dict-Salpetcr IMF and (4)ehosen IMF. Only NGC 2841 and NGC 2403 Seem to 
(5)-(6) for the Kroupa IMF. crj^g is the la-error on log{c2oo). M200 is iifoe galaxies whosc halo parameters are in full agreement 
^ ^ ^'S- with cuspy ACDM Einasto halos 



NGC925 

NGC2366 

NGC2403 

NGC2841 

NGC2903 

NGC3031 

NGC3198 

IC2574 

NGC3521 

NGC3621 

NGC4736 

DD0154 

NGC5055 

NGC6946 

NGC7793 



0.71 
1.01 
0.96 
1.04 
1.58 
0.88 
0.83 
0.70 
0.61 
0.78 
0.77 
0.89 
0.54 
0.75 
1.12 



46.9 
0.05 
0.06 
0.08 
0.04 
0.05 
0.04 
0.19 
1.73 
0.11 
1.35 
0.06 
0.06 
0.53 
0.01 



10.07 
9.49 
11.48 
12.44 
11.59 
10.24 
11.14 
10.66 
11.15 
11.59 
9.48 
9.97 
11.10 
10.72 
10.00 



0.78 
1.02 
1.00 
1.41 
1.72 
1.22 
0.95 
0.72 
1.04 
0.69 
1.40 
0.90 
0.72 
0.95 
1.18 



0.03 
0.05 
0.05 
0.02 
0.06 
0.02 
0.02 
0.14 
0.07 
0.13 
0.07 
0.06 
0.14 
0.02 
0.04 



10.15 
9.51 
11.50 
12.34 
11.61 
10.62 
11.16 
10.66 
10.99 
11.83 
9.92 
9.96 
11.19 
11.63 
10.18 



• E 



avarro 



et al.l 



ACDM dis sipationless simulations (e.g 
[200l[2010l) . 

A first caveat to this comparison is that, despite the 
wide global range of halo masses explored in these simu- 
lations, the halo masses do not span the whole range of 
our galaxies, namely 10^-5~12.5 A larger mass 

range of simulated galaxy-sized halos is required for fu- 
ture comparisons with observations. 

Secondly, and more importantly, these simulations usu- 
ally ignore the gaseous component of galaxies, whose dis- 
sipative nature leads them to concentrate in the cores of 
halos and possibly force the dark matter to respond to 
the concentration of the dominant central baryons. For 
this reason, we have also compared our results with those 
from lTissera et all (j2010( ) who have included the physics 
of bar yons in the galactic halos of the Aquarius simula- 
tions (Springel et al.ll2008[) by means of hydrodynamical 
models. 

Finally, the innermost regions of galactic halos are still 
hardly resolved by cosmological numeric al simulations 
(at sca les r < r2oo/1000), as noted in, e.g.. lNavarro et aO 

We distinguish here between fits performed with a 
free index and fits done with a fixed index in order to 
investigate whether the typical Einasto index seen in 
ACDM simulations matches the observations. We have 
thus performed fits with a fixed index n = 6.2, which 
is the average value of galaxy-sized halos presented in 
iNavarro et all (|2004[) . 

6.1. Models with a free index 

Figures [6] and [7] show that most of the observed and 
simulated galaxy-sized halos share the same characteris- 
tic scale densities and radii but not their indices. The 
indices of half of our observed halos are usually less than 
half the mean value of the simulated ones (i.e. n < 3). 
The Einasto indices for ba ryons-|-dark matter simulations 
from lTissera et al.l (|2010D differ even more from our ob- 
servations than do the pure dark matter halo indices. 
The addition of hydrodynamics in the ACDM simula- 
tions drags more dark matter towards the centre of ha- 
los, as likely caused by gas dissipational effects, thus 



Another interesting point concerns the results for the 
few dwarf and low surface brightness galaxies DDO 154 
and NGC 2366. Due to their dark matter dominance 
they are ideal candidates for a direct comparison with 
ACDM low mass halos (10^ h'^ Mq). They perfectly 
agree with simulated dwarf galaxy halos in the density- 
radius parameter space (the two closest symbols to the 
red triangles in Figs [6] and [7]) but they have a low Einasto 
index, contrary to all simulated dwarf-galaxy-sized halos. 

6.2. Cuspy Einasto models with a fixed index 

We now investigate whether Einasto halos with an in- 
dex close to the typical index seen in cosmological simu- 
lations agree better with the observations than the NFW 
and Iso halos. The Hi RCs have therefore been fit with 
a Einasto fixed value {n = 6.2) model. These fits have 
two free parameters, like the NFW and core halos. The 
results for both IMFs are reported in Tab. [5] and the RC 
decompositions are displayed in Fig. [12] (Kroupa IMF 
only). Figure [13] compares the quality of fits between 
the fixed Einasto index models to the free Einasto index 
models (left-hand panel), and to the NFW and Iso mod- 
els (right-hand panel), whereas Figure [H] displays the 
parameter space [p_2,r_2] obtained for n = 6.2. 

As seen in Fig. 1131 the fits performed with a fixed in- 
dex n = 6.2 are significantly worse than those with a free 
index, regardless of the IMF (84% and 89% of our sam- 
ple for the diet-Salpeter and Kroupa IMFs, respectively). 
They also provide worse results than the Iso model, re- 
gardless of the IMF (65% and 55% for the diet-Salpeter 
and Kroupa IMFs, respectively), and better results than 
the NFW model fits, regardfess of the IMF (75% and 80% 
for the diet-Salpeter and Kroupa IMFs, respectively). 

The parameter space of the characteristic scale densi- 
ties and radii still shows a correlation (Fig. [Hj). but with 
more scatter than seen in Figs. [6] and [7] for fits done with 
a free index. The striking feature is the larger range of 
densities and radii than for the fits performed with free 
indices. The relation between r^ and p(r_2) is steeper 
than the previous one given in H5.ll and Fig. [7] Besides 
the fact that extreme halos of low density and large ex- 
tent appear unrealistic for galaxies, the agreement with 
the expected galaxy-size range seen in cosmological sim- 
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Fig. 12. — Same as Fig. |4]but with a fixed index n = 6.2. 



ulations becomes worse and several fit halos exceed the 
cluster-sized regime. We therefore conclude that Einasto 
halos with a fixed index n ~ 6 are ruled out by the anal- 
ysis of RCs of galaxies. 

7. TWO-PARAMETER EINASTO MODELS 

In the previous section we have studied two-parameter, 
fixed index Einasto halos that were consistent with CDM, 
but noted they were not consistent with the observations. 
In this section, we try to identify other two-parameter 
Einasto halos and ask whether there are any Einasto 
models that can describe the observed halos better than 
the Iso or NFW halo models. Before answering these 
questions, we emphasize that our galaxy sample does not 
pretend to be representative of the whole diversity of spi- 
ral galaxies in the local Universe. Nevertheless, in this 
section, we aim at finding a family of generic Einasto ha- 



los that reproduce fairly well the RCs of THINGS spiral 
galaxies, at least as good as the usual core or cosmolog- 
ical cuspy models. 

7.1. Cuspy Einasto halos versus the NFW halo 

We estimate the critical index value where the NFW 
halo becomes better than the Einasto halo for the ma- 
jority of the current sample. We find n^p ^ 10 — 11 for 
the diet-Salpeter IMF and n^p ~ 11 — 12 for the Kroupa 
IMF. Those values have been derived by increasing pro- 
gressively the (fixed) index in the fits and by identifying 
the turnover index where more than half of the NFW 
fits from the sample becomes better than the Einasto 
fits. Note that those models represent halos which are 
less and less dense, and larger as the index increases. 
Another advantage of these models with respect to the 
NFW model is the better agreement with ACDM simu- 
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TABLE 5 

Fit parameters of the Einasto halo to THINGS galaxies for a fixed index n = 6.2 and mass-to-light 

RATIOS DERIVED USING A DIET-SALPETER (lEFT PANEL) AND KrOUPA (rIGHT PANEL) IMFS. 







diet-Salpeter 








Kroupa 






Galaxy 


Xr 


P-2 
(10-3 






r-2 
(kpc) 




P-2 

(10-3 pc-3) 




r-2 
(kpc) 


NGC925 


2.2 


1.2x10-5 ± 


8.0x10- 


7 


> 1000 


1.4 


1.8x10-'' ± 1.5x10" 




1.6xl05 ± 3.3xl05 


NGG2366 


0.6 


1.8x10-2 ± 


8.9x10" 




159.0 ± 80.9 


0.6 


2.0x10-2 ± 9.5x10" 




150.0 ± 73.8 


NGC2403d 


0.6 


5.8x10-1 ± 5.0x10- 




26.7 ± 1.5 


0.6 


1.2 ± 8.7x10-2 




17.8 ± 0.8 


NGC2403 


0.6 


5.5x10"^ ± 4.9x10" 




27.7 ± 1.6 


0.6 


1.1 ± 8.5x10-2 




18.2 ± 0.8 


NGC2841 


0.3 


3.6 ± 2.1x10^1 




18.9 ± 0.6 


0.3 


14.3 ± 7.7x10-1 




10.0 ± 0.3 


NGC2903 


0.3 


31.6 ± 2.5 






4.7 ± 0.2 


0.3 


48.1 ± 3.9 




3.9 ± 0.1 


NGG2976 


2.0 


4.7x10-5 ± 


1.3x10" 




> 1000 


0.5 


4.2x10-"' ± 3.5x10" 


4 


3.6x10^ ± 1.7x10^ 


NGG3031 


4.3 


1.2x10-2 ± 


2.2x10" 


2 


417.0 ± 760.0 


3.7 


14.4 ± 7.3 




5.6 ± 1.5 


NGC3198d 


1.9 


1.1x10-1 ± 


3.6x10" 


2 


62.9 ± 13.6 


1.8 


5.3x10-1 ± 1.3x10" 


1 


26.1 ± 3.6 


NGC3198 


3.2 


9.7x10-2 ± 


4.2x10" 


2 


66.9 ± 18.9 


2.6 


4.9x10-1 ± 1.5x10" 


1 


27.2 ± 4.7 


IC2574 


0.3 


2.4x10-* ± 


8.0x10" 


5 


> 1000 


0.3 


4.5x10-" ± 1.3x10" 


4 


9170 ± 4520 


NGG3521 


8.3 


4.2x10-^ ± 


1.5x10" 


6 


> 1000 


5.7 


2.6x10-1 ± 2.7x10" 


1 


50.3 ± 36.7 


NGG3621 


0.8 


3.9x10-2 ± 


6.7x10" 


3 


147.0 ± 19.4 


0.5 


3.2x10-1 ± 3.1x10" 


2 


37.6 ± 2.3 


NGC4736 


1.6 


3.3x10-1 ± 


1.1 




12.1 ± 24.8 


1.6 


78.8 ± 52.2 




1.3 ± 0.4 


DD0154 


0.8 


4.4x10-2 ± 


1.2x10" 


2 


54.1 ± 11.8 


0.8 


4.6x10-2 ± 1.2x10" 


2 


52.3 ± 11.4 


NGC5055 


8.3 


3.7x10-** ± 


3.2x10" 


7 


> 1000 


1.5 


6.7x10-3 ± 2.6x10" 


3 


541.0 ± 182.0 


NGC6946 


1.5 


6.2x10-5 ± 


5.0x10" 


5 


> 1000 


1.0 


1.1x10-1 ± 2.9x10" 


2 


106.0 ± 21.8 


NGG7331 


3.0 


6.4x10-^ ± 


6.7x10" 


7 


> 1000 


0.3 


9.2x10-2 ± 2.7x10" 


2 


112.0 ± 24.4 


NGG7793 


4.1 


1.1x10-1 ± 


5.3x10" 


2 


89.2 ± 39.1 


3.9 


2.9x10-1 ± 1.1x10" 


1 


42.8 ± 13.2 


NGC7793S 


1.5 


5.4x10-3 ± 


3.0x10" 


3 


> 1000 


1.8 


4.1x10-2 ± 1.8x10" 


2 


324.0 ± 162.0 
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Fk;. 13. — Comparison of reduced X'^ for Einasto halo models 
done at fixed index (n = 6.2) with Einasto halo models done at free 
index (left). Open and filled symbols are for the diet-Salpcter and 
Kroupa IMFs, respectively. Comparison of reduced for ri = 6.2 
Einasto halo models with the NFW and Iso halos (respectively 
filled and open symbols) with the Kroupa IMF (right). See Tab.|5] 
for results obtained with the diet-Salpeter IMF. 



lations. 

We repeat the same procedure but towards smaller 
indices and find niow ^ 0.5 (diet-Salpeter IMF) and 
niovj ^ 1 (Kroupa IMF) as the values where the ma- 
jority of constrained Einasto fits turn better than NFW 
ones. This result is completely expected because small 
indices provide core halos, whose halo shape is naturally 
better for the sample than any cusp. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
log(r_2 [KPC]) 

Fig. 14. — Parameter space of Einasto halos with a fixed index 
(n = 6.2) for the Kroupa IMF. See Tab. [5] for results obtained 
with the diet-Salpeter IMF. Results for the few halos having an 
extremely large scale radius are not shown for clarity. Colored sym- 
bols correspond to the same cosmological halos as in Figs. Eland [7] 
Additional blue, open triangles corre spond to simulated cluster- 
sized halos from [Navarro et al.l II2004I) . A dashed line is a power 
law fit to the relationship log(p_2) = (—2.25 ± 0.04) log(r_2) -|- 
(0.16 it 0.05). As a comparison, a dotted line represents the linear 
fit to the relations hip o btained from Einasto models done at free 
Einasto index (see 115. H and Fig. [3. 



7.2. Core Einasto halos versus the pseudo-isothermal 

sphere 

In contrast with the cuspy halos, we find the cored 
(n < 4) halos at the low end of the Einasto index dis- 
tribution. We investigate the goodness of fit of two- 
parameter Einasto models performed with a fixed, small 
index to find a critical index where the majority of such 
constrained Einasto fits become worse than those per- 
formed with the Iso model. It turns out a low limit can- 
not be obtained for the diet-Salpeter IMF; the Iso model 



is always a better fit than the constrained Einasto model 
at a level of ^-^70%, which is likely explained by the rel- 
atively scattered index distribution, as seen in i)5.1l In 
other words it is almost impossible to find an index that 
would satisfy most of all halos at the same time for that 
IMF. It is not the case for the Kroupa IMF for which 
"■low ~ 2 indicates the turnover index. Indeed, 65%, 55% 
and 45% of all the sample is better fit by the Iso model 
than n = 1, n = 2 and n = 3 (respectively) Einasto fits. 
Notice this value appears consistent with the mean index 
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n = 1.4 ± 0.8 of non-euspy halos fit under the Kroupa 
IMF assumption (^EUl- 

For larger indices n^p '--^ 5 is the limit where most 
of constrained Einasto fits become worse than Iso ones. 
Past this value, the core nature of the Iso model starts to 
dominate the cuspy nature of these Einasto halos. The 
index range 2 < n < 4 of constrained Einasto core mod- 
els is therefore small compared with that of constrained 
Einasto cusps. 

7.3. Towards a new two-parameter density profile 

Both results of i jT.ll and ij7.2l make us conclude that the 
family of two-parameter Einasto models with 4 < ti < 12 
is more suited for the modeling of observational cuspy ha- 
los than is the NEW model, whereas the Einasto model 
with 2 < n < 4 is more suited for the modeling of cored 
halos than is the isothermal model. While allowing the 
variation of the Einasto index is essential in finding the 
best fit models to the THINGS RCs (gl gS]and P, those 
fits at fixed indices point out the exponential shape of the 
density profile also contributes in better modeling the 
mass distribution of the galaxies. We thus suggest that 
a new shape of dark matter density, where the volumic 
density profile decreases exponentially as a function of 
radius following Eq. [31 may be a good alternative to the 
usual mass divergent, pseudo-isothermal and NEW mod- 
els, with Einasto indices chosen to match the particular 
type of halo (core or cusp) one wants to model. 

8. CONCLUSIONS 

Eor the first time, the Einasto dark matter halo has 
been applied to model the mass distribution of nearby 
and undisturbed galactic disks. Our results show that 
the Einasto model is a good fit of the dark matter con- 
tribution to highly-resolved Hi rotation curves of galax- 
ies, better than other models like the (cored) Iso model 
or the (cuspy) NEW model. The significant improve- 
ment is caused by the parameter that shapes the radial 
behaviour of the mass density profile, the Einasto in- 
dex. The big advantage of the Einasto formalism over 
the NEW cusp or the pseudo-isothermal sphere lies on 
its third parameter, the index, which allows for extra 
freedom in describing different shapes of density profiles. 
The larger (respectively, smaller) is the Einasto index, 
the cuspier (shallower) is the halo in its central parts. 
Since there is no unique index describing all halos, it is 
not possible to scale the density profile from one galaxy 
to another with a different mass. In other words, we find 
no evidence for a universal dark matter Einasto halo, at 
least in the galaxy sample we have analyzed. This re- 
sult corroborates the one found in dissipationless ACDM 
simulations. Noteworthy however is the fact that the in- 
dex appears to scale with the halo mass for halos more 
massive than M200 ^ 2 x 10^^ Mq. 

We find that the Einasto index is generally smaller 
than the one inferred from cosmological simulations, by 
a factor of two or more. This discrepancy between ob- 
served and simulated Einasto indices extends the core- 
cusp controversy (Iso versus NEW) to the Einasto for- 



malism. A ACDM Einasto halo is cuspier than an ob- 
served galactic Einasto halo. In our sample, galaxies 
with a cuspy (i.e. with n > 4) Einasto halo are not fre- 
quent. It is interesting to note that our sample consists 
of both normal, bright galaxies and of a few dwarf, low 
surface density galaxies, while the cusp-core controversy 
usually deals with small and low surface density galaxies 
only. Cored dark matter halos are thus found in galaxies 
whose baryonic content significantly differs, for halos less 
massive than M200 2 x 10^^ Mq- 

The RC decomposition has also shown that cuspy 
Einasto halos with fixed index 4 < n < 12 give a better 
description of RCs than does the NEW model. Con- 
versely core Einasto models with fixed index 2 < n < 4 
give better RC fits than the Iso model. It implies that the 
exponential decrease of the Einasto dark matter density 
is also responsible for the improvement in modeling the 
mass distribution of galaxies, in addition to the effect 
of the Einasto index. Those constrained Einasto mod- 
els can be a good alternative to Iso or NEW models for 
galactic dynamics, but we re-iterate here that Einasto 
models with free index are preferable over any of those 
two-parameter models. 

The a nalysis has shown that the stellar masses derived 
with the lBell fc de Jon3 ()2001f ) prescription under a dict- 
Salpeter IME hypothesis are in conflict with the Einasto 
model. This is not the case for the Kroupa IME. In fu- 
ture papers from this series, we will perform mass models 
with free stellar masses and with other fixed stellar mass 
scalings, like the one described in iChabrier (20031). Addi- 
tional improvements of dynamical models will take into 
account the molecular gas content, where available. 

Eurthermore we will investigate possible correlations 
between the properties of Einasto dark matter halos and 
baryons. The mass distribution for a large sample of 
low surface density disks will also be studied. We expect 
that those galaxies exhibit low Einasto indices because 
of the shallow nature of the density profile in their core, 
but to which extent their values compare with those of 
normal, higher surface brightness disks like most of our 
current THINGS subsample is not known yet. Both these 
samples will help us to measure the variation of dark halo 
inner shape as a function of galaxy mass. 

We finally emphasize the importance of testing the 
Einasto halo model for objects at extreme ends of the 
galaxy luminosity function (dwarf spheroidals and giant 
ellipticals) . This will be crucial for getting a complete de- 
scription of the observational behaviour of Einasto mod- 
els over several decades of halo mass and comparing with 
expectations from ACDM simulations at similar cosmo- 
logical scales. 
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